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MULTIPATH  IN  THE  THREE-DIMENSIONAL  UNDERWATER  ARRAY 

A  planar  array  of  widely  dispersed  hydrophones  deployt’d  in  a  horizontal 
plane  several  hundred  meters  below  the  surface  of  the  sea  been  analyzed, 
and  reported  earlier  [1,  2].  Because  of  the  dispersive  nature  of  the  under¬ 
water  medium, signal  energy  approaches  the  array  from  a  source  along  a  number  of 
refracting  paths.  This  resul ts  in  multipath  interference  a^ci  a  possible 
consequent  loss  of  output.  Ray  arrivals  from  the  source  typically  fall  into 
a  range  of  vertical  angles  at  the  receiver  which  are  t  10°  relative  to  the 
horizontal.  The  vertical  beamwidth  of  a  planar  antenna  is  large  enough  to 
accept  all  rays  in  such  a  range,  hence  it  is  multipath  sensitive.  A  three- 
dimensional  array  is  capable  of  a  sharper  vertical  focus  and  will  be  lens  sensi¬ 
tive  to  the  kind  of  multipath  typical  in  this  application.  in  fact,  by  suitably 
processing  the  array  output  there  is  the  possibility  that  the  multipath  arrivals 
can  be  separately  received  and  then  combined  in  phase  to  achieve  an  "angle  cf 
arrival  diversity1*  system, ns  suggested  in  Figure  3.1.  Such  systems  have  been 
proposed  for  tropospheric  sl alter  receivers.  . 

We  analyze  the  mean  properties  of  such  a  system  below.  If  the  array  is 
focused  to  look  in.  the  y-z  plane  its  respond  to  a  signal  at  angle  y  relative 
to  the  plain.*  is 
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FIGURE  3.1  MULTIPLE  VERTICAL  BEAMS  FOR  MULTIPATH  RESOLUTION 
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Here  G  is  the  vertical  angle  to  which  the  array  is  to  be  focused; 


(x  ,  y  ,  2  )  is  the  position  of  the  nth  array  element.  (1)  assumes  M  ray 
n  n  n  * 


arrivals  each  given  by  B  e  with  vertical  arrival  angle  0  . 

m  it. 


We  concentrate  on  tiie  output  when  the  source  is  on  the  main  beam;  that 
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We  calculate  the  mean  power  response  of  the  array  given  by 


M  M  N  N 


/■,  s  l  1  11  «  V 

<|A  (j,0s)|^  -  l  l  l  l  <».  B  > 


>  Ill-  ill.,' 

ml  m2  nl  n2  1  * 


jkfy  (sin0  -  sinO  )  -  y  (sin9  -  sinO  )] 
n, 


V 


s  n2  m2  S 


jkfz  (cos0  -  cos9  )  -  z  (cos9  -  cos0  )] 

s  n2  ^  s  ^ 


j  (4>  -  ) 

/  ml  m2  \ 
*  <£  / 


(3) 


We  have 


j  (v>  -  o  ) 

m  m2  v 

<6  / 


mln2 


(A) 


the  Kronecker  delta,  so  that 


<ia2(|,v:>  -  ?  I.  !<»;><c  '• 

m=l  n^i  n  =1  N 


jk(y  -y  )(sinO  -sinO  ) 
n  1  n  2  m  s 


ik(z  -?  )  (cosO  -cosO  ) 

m  s 


> 


M  .  2.  . 


m  m  jk  y  (sinO  -  sinO  ) 
N  N  ;n,  N  m  s 


l<V>|»*  l  £<•  1 

m=l  '  l_  n  =1  n„  =  l  x 


n]^  n2  jk  z  (cos'1  -  cosO  ) 

m  s  v 

*  e  / 


•  A 


-jkv  (sinO  -  sinO  ) 


“ikz  (cosn  -  cosO  ) 
n  2  m  s 


(5) 


oru  No.  n 


We  assume  all  vectors  (v  ,  z  ),  n.  *  1,  2,  ...N  identically  distributed. 

'  n .  n .  i 
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Furthermore,  the  random  variables  y  and  2  are  assumed  independent  and 
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symmetrical  around  the  origin.  Then 
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The  expectations  oa  the  right  are  characteristic  functions, 


$(jt)  =  <%jtu  (7) 

where  the  random  variable  u  is  either  v  or  z  and  t  is  correspondingly 
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either  k(sinG  -  sinO  )  or  k(cosO  -  cosO  ).  The  random  variables  will  here 
ms  ms 

be  specified  as,  cither,  uniformly  distributed  in  an  interval  (~h,  h) ,  or 
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normally  distributed  around  zero  with  variance  0  .  Thus  for  the  uniform 
case 
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and  for  the  normal  case 
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Jf  the  variables  y  and  z  arc  both  normal  with  variance  0,,  and  0,, 
J  n  n  y  2 

respectively,  we  have 
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If  the  variable  y  is  normal  with  variance  0  and  the  variable  z  is 
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uniform  in  (-h,  h)  then 
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An  inspection  of  (10)  or  (11)  leads  to  the  conclusion  that  if  the 
vertical  dimension  of  the  array  is  in  the  order  of  10  wavelengths  the 
vertical  bear.width  will  be  about  t  1°.  Flirt  lerr.ioro  rays  entering  through 
this  narrow  beamwidth  (hence  excluding  c  t  h  e  r  rays  arriving  at  vertical 
angles  outside  the  vertical  beantwidth)  will  no  sufficiently  compact  to 
avoid  the  effect  of  phase  decorrelation  across  the  array. 

We  are  therefore  led  to  propose  the  following  concept.  Let  the 
array  simultaneously  form  contiguous  vertical  beans  as  indicated  in 
Figure  3.1.  Outputs  corresponding  to  each  beam  will  be  simultaneously 
present.  Ihese  outputs  are  then  coherently  combined.  The  operations 
required  are  as  indicated  in  Figure  3.2.  The  mechanism  being  suggested 
is  similar  to  that  used  in  angle  of  arrival  diversity  communication  sys¬ 
tems  with  maximal  ratio  combining  of  the  diversity  signals.  As  a  rule 
in  these  systems  each  diversity  branch  has  a  separate  directive  sensor 
and  preamplifier.  Here  sensors  and  preamplifiers  are  common  for  all 
b  ranchos . 
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Because  sensors  are  common  one  may  question  the  effect  of  noise 
generated  at  the  sensor  or  preamplifier  input.  Will  such  noise  be 
independent  when  observed  at  the  point  of  combination  of  the  diversity 
branches?  The  following  is  a  discussion  of  that  peint. 

A  filter 


H (f )  -  A(f )  e^(f) 

which  acts  as  a  constant  gain  device  and  constant  phase  shifter  -  that  is, 
with 


A(f)  =  A 

<J>(f)  -  -v  sgn  f  ($  a  constant) 

can  be  represented  by 

H(f)  =  A  e  ^?S",n  ^  =  A  ces  (-:  sgn  f)  +  j  A  sin(-0  sgn  f) 

=  A  cos:  -j  A  sin  ;  sgn  f 

A  filter  with  frequency  characteristic 

H  (f)  =  “j  sgn  f 

is  a  Hilbert  transforming  filter  so  that  a  wave  function  n(t)  applied  to 
H(f)  as  defined  above  emerges  as 

n  (t)  -  A  cos,*'  n(t)  +  A  sin*  n(t) 
o 

where  n(t)  is  the  Hilbert  Transform  of  n(t). 

By  direct  application  of  the  definitions  and  by  use  of  the  statistical 
properties  of  the  Hilbert  Transform  one  can  show  that  for  stationary, 
zero  mean,  processes, 
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^ (t+i)  nQ(t)^>  =  A2<^n(t+i)  n(t)^>  =  A2  Rr(t) 

Input  and  output  autocorrelation  functions  are  proportional.  R^(t) 
the  input  autocorrelation  function.  Also, 

/  n  (t+t)  n(r)  N  -  A  cos <t>  R  (t)  +  A  sin<}>  R  (t)  , 
x  o  /  n  n 

giving  a  relationship  between  input-output  cross-correlation  function  and 

the  input  autocorrelation  function.  R  (t)  is  the  Hilbert  Transform  of 

n 

R^  ( t) .  Finally  for  an  input  n(t)  applied  to  separate  filters 
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the  cross-corrclat i an  function  of  the  two  outputs  n  (t)  and  n  (t)  is 
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Consider  now  the  block  diagram  of  Figure  3.2.  Each  branch  is  com¬ 
prised  of  the  sum  of  N  inputs,  one  from  each  array  element  and  phase 
shifter.  Branch  1,  for  instance,  contains  a  wave  function 


N  (t)  =  l  n  ,(t)  =  l  IcosJ  n  (t)  +  sin*  n  (O] 
1  n=l  onl  n=l  nl  n  nl  n 


where  n  (t)  is  the  output  at  the  n’th  sensor,  and  n  (t)  is  the  phase- 
n  onl 


shifted  output  of  the  n’th  sensor  contributing  to  branch  1.  .-*>  is  defined 

Tnl 

in  Fig, ii re  3.2. 


QFR  No.  31 


The  output  N^(t)  will  ultimately  be  applied  through  the  second  filter 

-jtf'jSgnf 

with  charac teristic  given  by  A^e  resulting  in  an  output 
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For  all  I  branches  taken  together  we  get 
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We  view  the  n  (t)  to  bo  Gaussian  noise,  present  at  the  sensor  outputs 
n 

generated  in  the  sensor,  its  preamplifier  and  trie  sensor's  immediate 

surroundings.  We  assume  n  (t)  is  independent  of  n  (t)  for  n  ^  in.  Thus 

n  fit 

we  exclude  external  noise  which  may  !>V’  correlated  across  several  sensors 
The  variance  of  N(t)  is  then 
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of  n  (t)  and  n  (t)  for  all  n  and  m  for  a  Gaussian  process.  It  can  be 
n  m 

shown  that 


<nn2U)>  -  <s>> 


so  that 


<N2(t))>  =  Jx<nn2(t)>  <A.A.  cosWn.-  4>nj  +  *.)> 


To  continue  this  analysis  ve  require  the  joint  statistical  properties 
of  the  phase  shifts  and  the  amplitude  factors.  For  our  purposes  at  present 
we  may  assume  the  constant  for  all  i  »  1,  2,  ...1.  But  the  phase  shift 
properties  are  needed,  Note  the  following.  If  =  1,  all  i,  and 

<  cosKr  'V  h  '  VV  °* 

except  when  i  =  j ,  then 

<N2(t)>  -  I  l  <ry2(t>> 

However,  if  the  angular  differences  were  small  so  that 


=  1 


for  all  j  and  j  then 

<k2<o  }  ■  l2l  <"/<'>> 

In  t no  latter  case  the  branch  noises  are  correlated  and  add  coherent lv. 
In  the  furrier  case  they  are  uncovrelated  and  add  incoherently . 

The  difference  angle  statistics  are  under  investigation  and  will  be 
reported  later.  A  preliminary  calculation  indicates  that  for  the  array 
size  envisioned  in  this  application  l lie  angular  djfft’rences  may  be  large 
enough  feu  the  tirst  rendition  above  to  he  approximately  correct. 
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